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This chapter starts with a summary of the atomistic processes that occur 
during epitaxy. We then introduce density functional theory (DFT) and 
describe its implementation into state-of-the-art computations of complex 
processes in condensed matter physics and materials science. In particular 
we discuss how DFT can be used to calculate parameters of microscopic 
processes such as adsorption and surface diffusion, and how they can be 
used to study the macroscopic time and length scales of realistic growth 
conditions. This meso- and macroscopic regime is described by the ab initio 
kinetic Monte Carlo approach. We discuss several specific theoretical stud- 
ies that highlight the importance of the different diffusion mechanisms at 
step edges, the role of surfactants, and the influence of surface stress. The 
presented results are for specific materials (namely silver and aluminum), 
but they are explained in simple physical pictures suggesting that they also 
hold for other systems. 
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1. INTRODUCTION 

Many basic concepts for surface diffusion and crystal growth have been 
developed already more than 40 years ago (see e.g. |I]-£3j). In recent 
years the subject has attracted significant attention that is largely due 
to new experimental advancements, progress in theoretical methods, and 
the importance of crystal growth for technological applications. A thor- 
ough knowledge of the motion of atoms at surfaces is a key factor to the 
understanding of chemical reactions at surfaces, of crystal growth, and to 
the question under which conditions thermal equilibrium can be achieved 
(at least locally) at crystal surfaces. In typical growth experiments de- 
position rates are of the order of a monolayer (ML) per minute, and the 
diffusion of atoms at the surface is too slow (at least for some processes) 
so that thermal equilibrium is often not reached. As a consequence, the 
structures which occur at surfaces are usually a result of the kinetics. 
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Figure 1. Three different growth modes of heteroepitaxial growth. 



Nevertheless, under certain conditions the resulting structures are ruled 
by thermal equilibrium, i.e., they correspond to the minimum of the free 
energy. This may occur when the adatom mobility is very high, the depo- 
sition rate is low, or when growth is interrupted and the sample annealed. 
Concerning the acquired surface morphology one then distinguishes the 
following three different growth modes (cf. Fig. [I]): Frank- van der Merwe 
or layer- by- layer growth, Volmer- Weber growth, where three-dimensional 
islands are formed and the overlayer does not completely cover the ex- 
posed substrate surface, and Stranski-Krastanov growth with layer-by-layer 
growth supplanted by island growth. Bauer || had discussed the conditions 
for these growth modes and pointed out that the actually realized surface 
morphology (under thermal equilibrium conditions) is controlled by the 
competition of the surface energy of the substrate, the surface energy of 
the film, and the interface energy of the film and substrate. 

The modern treatment of growth under non-equilibrium conditions starts 
with the seminal paper by Burton, Cabrera, and Frank (BCF) [f3j, who 
realized that a surface usually is not perfectly flat but has imperfections 



such as steps. Steps are the boundary between regions that correspond to 
upper terraces and regions that correspond to lower terraces. They might 
occur at random or can be created in a controlled manner by cutting the 
surface in an orientation close to a low index plane. In the latter case the 
surface is called a vicinal surface. At steps and particularly at kink sites 
adatoms are bound most favorable and BCF assumed that on a stepped 
surface growth occurs by attachment of deposited adatoms to steps that 
subsequently advance. This growth mode is called step flow and is indeed 
easily achieved experimentally at higher temperatures. It results in flat 
films. The step-flow growth mode requires a high mobility of the deposited 
atoms, so that they reach the existing steps before meeting other adatoms. 
Thus this situation is close to thermal equilibrium. 

Usually it is desirable to grow crystals at not too high temperatures. 
Then the adatom mobility is lower, and deposited atoms might not reach 
an intrinsic step edge. Instead, they will wander around on the surface and 
meet other atoms, eventually nucleating a new small island. Further atoms 
will be caught at the edges of these islands (or create new ones), so that 
growth proceeds via two-dimensional growth of the islands. Additionally, 
during the initial stage of film growth the chance that newly deposited 
atoms land on an island is small. Nucleation and growth of islands can 
be described by phenomenological rate equations [|5],£| and we will discuss 
this approach briefly in the following Section. The competition between 
step flow growth and growth via nucleation, spread, and eventual coales- 
cence of islands on the terraces of a vicinal surface can be captured if one 
incorporates the rate equation formalism into the BCF equations [jTj . 

Using rate equations one usually describes diffusion by effective param- 
eters but unfortunately lacks a detailed understanding of the microscopic 
mechanism behind them. For example, it has been discussed by several au- 
thors [|8|— [T^jl that surface diffusion can occur via two different mechanisms: 



i) An adatom may simply hop from one low-energy site to another one, 
while the substrate reacts only modestly by local relaxations, or 

ii) an adatom may diffuse by atomic exchange where it changes place 
with a substrate atom and the ejected substrate atom moves further. 

These two mechanisms are operative at the flat regions of the surface, but 
also for diffusion across steps or parallel to step edges. The interplay of 
those different diffusion processes significantly affects the shape of growing 
islands. We will also see an interesting dependence of these two mechanisms 
on surface stress, that exists at free surfaces and also results from lattice 



mismatch during heteroepitaxy. 

Kinetic limitations might lead to either two-dimensional growth or three- 
dimensional growth. The former is achieved when atoms that land on 
islands can easily move down. However, when atoms which land on an 
island are hindered to move down, islands nucleate on top of islands, and 
a three-dimensional structure results (see also Ref. [|T5fl ). As in the above 



discussion of thermal equilibrium one also often labels the kinetic growth 
modes according to the surface morphologies of Fig. [T], although the kinetic 
and the thermodynamic limit represent totally different physics. 

In order to gain insight into growth phenomena it is necessary to ex- 
amine all possibly relevant microscopic processes on the atomic scale. In 
Section [27T] these processes will be identified. Section [Z72] then summarizes 



some aspects of the description of growth by rate equations, and Section |2T3| 
analyses the conditions at which two-dimensional growth is attained. 

Theoretical methods that are promising for a reliable description of sur- 
face diffusion and growth are described in Section [3]. In particular we give a 
brief review of density functional theory in Section \5.2\ . Then, in Section |3.3| 
we sketch how DFT is implemented in actual computational schemes and 
utilized to identify microscopic processes and to obtain growth-relevant 
parameters. Section £I7t| describes how these ab initio parameters can be 
used to predict or analyze the temporal and spatial evolution of epitaxial 
growth on macroscopic scales. Ab initio kinetic Monte Carlo simulations 
make the connection between the atomic scale and time and length scales 
of realistic growth conditions. 

Section f| then presents some recent results of close-packed fee surfaces. 
We start with an analysis (but all these results are still predictions) of the 
aluminum and silver (111) surfaces (Sections ffTT] and 4.2). Sections [O] 
and |4.4| discuss results for aluminum and silver (100). In particular, we 



consider the effects of stress for the two silver surfaces (Sections \i:z.\\ and 
4.4.2[ ), possibly modified by strain due to heteroepitaxy. For the Ag(lll) 



surface we discuss in Section |4.2.2| how adatom motion, the island density, 
and consequently the growth mode can be influenced by surface active 
contaminants, so-called surfactants. 

2. ATOMISTIC PROCESSES AND RATE EQUATIONS 

2.1. Atomistic processes 

The conceptually simplest growth technique is molecular beam epitaxy 
where an atom that has landed at the surface may either stay on the 
surface and wander around, or evaporate back into the gas phase. The 



latter happens at a rate T ev = rQ V exp(— E^/k^T), where T™ is the effec- 
tive attempt frequency (in some works this term is approximated by the 
vibrational frequency of the isolated adatom on the surface), E &( \ is the 
adsorption energy of the adatom, k& the Boltzmann constant, and T the 
substrate temperature (for simplicity we assume here a situation of atomic 
desorption, i.e., no formation of molecules). Typically E^ is larger than 
the activation barriers for other processes that occur on a surface so that 
regardless of the exact magnitude of Tq v (typically of the order of 10 13 s _1 ) 
evaporation can be neglected during growth. 

The different atomistic processes encountered by adatoms are illustrated 
in Fig. g. After deposition (a) atoms can diffuse across the surface (b) and 
will eventually meet another adatom to form a small nucleus (c) or get 
captured by an already existing island or a step edge (d). Once an adatom 
has been captured by an island, it may either break away from the island 
(reversible aggregation) (e) or remain bonded to the island (irreversible 
aggregation). An atom that is bonded to an island may diffuse along its 
edge (/) until it finds a favorable site. As long as the coverage of adsorbed 
material is low (say 6 < 10 %), deposition on top of islands is insignificant 
and nucleation of islands on top of existing islands practically does not 
occur. However, if the step down motion (g) is hindered by an additional 
energy barrier, nucleation of island on top of islands becomes likely (h). 

In principle it is possible that not just single adatoms but also dimers 
and bigger islands migrate (i). For example, a dimer might diffuse by the 
two atoms rotating around each other. Moreover, compared to a single 
adatom, a dimer may be less bounded to the substrate since the electrons 
of the two adatoms participate to the adatom-adatom bond and not only 
to the adatom-substrate bonds. Therefore, it may be expected a low ac- 
tivation barrier for the diffusion of dimers, but there is no clear evidence 
yet available. Finally, it is sometimes believed that a large island is com- 
pletely immobile. However, results of Wen et al. [|TBj for Ag/Ag (100) show 



that even large scale clusters with 10 2 to 10 3 atoms can diffuse at room 
temperature. Diffusion of a cluster can either happen by consecutive edge 
diffusion of single atoms from one side of the cluster to the other, or by 
some concerted motion of all atoms in the cluster. The importance of the 
diffusion of dimers or large islands during growth is an issue that deserves 
more attention in future research but will not be addressed any further in 
this chapter. 

Processes such as attachment to and detachment from step edges depend 
quite sensitively on the local environment, because chemical bonding is a 
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Figure 2. The different atomistic processes for adatoms on a surface: (a) deposition, 
(6) diffusion at fiat regions, (c) nucleation of an island, (d) diffusion towards and capture 
by a step edge, (e) detachment from an island, (/) diffusion parallel to a step edge, (g) 
diffusion down from an upper to a lower terrace, (h) nucleation of an island on top of an 
already existing island, and (i) diffusion of a dimer (or a bigger island). For the processes 
(a), (c), (g) and (h) also the reverse direction is possible, but typically less likely. 



rather local phenomenon and largely determined by the coordination of 
nearest neighbor atoms. Figure |3| displays three important geometries of 
step-edge atoms. In a bond cutting model of metallic bonding the energy of 
an atom scales as the square root of the local coordination (see Section |3TT| 
for more details). In this approach it follows that the binding energy of 
an atom at a kink site (atom 1 in Fig. £3]) equals the cohesive energy. This 
is indeed a general results and implies that kink sites help to establish 
thermal equilibrium of the surface with the bulk (note that a kink atom 
which detaches from a step creates a new kink in the step edge). Compared 
to an isolated adatom on the surface, the binding energy of an atom in the 
step edge (atom 2 in Fig. 0) is about 30 % and 60 % larger than the binding 
energy of atom 1 and atom 3, respectively. The above discussion refers to 
metallic bonding and is not necessarily valid for systems which can form 
covalent bonds and therefore prefer a certain low coordination. 

At not too high temperatures atoms will usually not detach from an 
island but diffuse along the edge. Eventually they will reach a higher 
coordinated site such as kink site 1. In general, fast edge diffusion leads 
to rather compact island shapes, but when the edge diffusion is strongly 
hindered, adatoms remain at the edge site where they reach the island and 




Figure 3. Three important geometries of atoms at step edges at an island on a fee (100) 
surface. 



islands acquire a fractal or dendritic form ("hit and stick mechanism" [|T7j] ). 



At a fee (111) there are two close-packed steps (see Fig. |4j), and be- 
cause more open steps typically have a higher step formation energy, these 
close-packed steps are expected to dominate the periphery of islands. The 
steps are labeled as {100} and {111} facets, referring to the plane passing 
through the atoms of the step and the atom of the substrate (often these 
steps are labeled A and B). Because of the microscopic difference of the two 
types of steps the diffusion along them will be different as well. It has been 
observed for growth of Pt on Pt (111) |TEj] that at a certain temperature 



the shape of the islands observed is triangular. At a higher temperature, 
triangular islands are observed again, but the triangles are rotated by 60°. 
More precisely, the islands are bounded only by {100}-faceted steps at a 
lower temperature while at a higher temperature the islands are bounded 
by {lll}-faceted steps. It was proposed |I8|] that this is a consequence 
of the different diffusion constants for migration along the two steps and 
particularly their different temperature dependences. 

The key idea behind the kinetic description of the growth phenomena is 
that processes occurring during growth, such as diffusion or desorption, are 
described by rates. The rate of a microscopic process j that occurs during 
growth usually has the form |jl3|-plj| 

T^ = —exp(-AF^/k B T) , (1) 

h 

where AF^' is the difference in the Helmholtz free energy between the 
maximum (saddle point) and the minimum (equilibrium site) of the poten- 
tial curve along the reaction path of the process j. T is the temperature, 
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Figure 4. The two different types of close-packed steps on a fee (111) surface. 



/cb the Boltzmann constant, and h the Planck constant. The free energy 
of activation AF^' needed by the system to move from the initial position 
to the saddle point is given by 



AF {j) = E, 



M 



U) 



TAS^ 



(2) 



Here E^ is the sum of the differences in the static total and vibrational 
energy of the system with the particle at the minimum and at the saddle 
point, and AS^l is the analogous difference in the vibrational entropy. The 



rate of the process j can be cast as follows: 
r«> = lf» exp(-E l />/k B T) , 



(3) 



where Iq = (k^T/h) exp(AS'|jy &b) is the effective attempt frequency. 
In the case of isotropic motion of an adatom on the surface it follows 
from Eq. (§) that the diffusion constant is D = D exp(-E%> /k B T) j22j. 

The prefactor Dq = l/{2a)T(j I 2 where / is the jump length and a the 
dimensionality of the motion (a = 2 for the surface). 



The two basic quantities in Eq. (§) are the attempt frequency Iq and 

the activation energy E^. Transition state theory (TST) [|20| , |ZT| allows an 

(j) 
evaluation of Tq within the harmonic approximation: 

d _ M U J (4 ) 

where Vj and z/* are the normal mode frequencies of the system with the 
adatom at the equilibrium site and at the saddle point, respectively, and 
3iV is the number of degrees of freedom. The denominator in Eq. (|j) 
contains the product of only 3iV — 1 normal frequencies, because for the 
adatom at the saddle point one of the mode describes the motion of the 
particle toward the final site and has an imaginary frequency. TST is only 
valid when E& is larger than k^T. 

The attempt frequency Tq shows a much weaker temperature dependence 
than the exponential and for typical growth temperatures it is of the order 
10 12 — 10 13 s -1 . When the barriers for two different diffusion events are 
different a compensation effect |23j may occur, i.e., Tq is larger for processes 



with a higher energy barrier. Indeed, a higher energy barrier usually implies 
a larger curvature of the potential well around the equilibrium site of the 
adatom. The corresponding vibrational frequencies of the adatom in such 
a potential are larger as well, which implies [see Eq. (||)] that the attempt 
frequency increases. 

To define and determine i? d (and other quantities important for the 
description of growth such as adsorption energies) we need to calculate 
the ground-state total energy of the adsorbate system for a dense mesh of 
adatom positions. This yields the so-called potential-energy surface (PES) 
which is the potential energy experienced by the diffusing adatom, 

£ PES (X ad , y ad ) = mm £ tot (X ad , y ad , Z ad , {Rj}) , (5) 

where i? tot (X ad , y ad , Z ad , {R/}) is the ground-state energy of the many- 
electron system (also referred as the total energy) at the atomic configura- 
tion (X ad , y ad , Z ad , {R/}). According to Eq. (g) the PES is the minimum of 
the total energy with respect to the z-coordinate of the adatom Z ad and all 
coordinates of the substrate atoms {R/}. Assuming that vibrational effects 
can be neglected, the minima of the PES represent stable and metastable 
sites of the adatom. Note that this PES refers to slow motion of nuclei 
and assumes that for any atomic configuration the electrons are in their 
respective ground state. Thus, it is assumed that the dynamics of the 



10 



initial transition final 

geometry state geometry 



.a 




(b) 





110] 



112] 



Figure 5. Diffusion via hopping (a) and exchange (b) on a fee (100) surface and diffusion 
along a {111} step on a fee (111) surface via exchange (c). 



electrons and of the nuclei are decoupled. This is the Born-Oppenheimer 
approximation that for not too high temperatures is usually well justified. 

Now consider all possible paths / to get from one stable or metastable 
adsorption site, R a a, to an adjacent one, R a d'- The energy difference E& 
between the energy at the saddle point along / and the energy of the initial 
geometry is the barrier for this particular path. If the vibrational energy is 
negligible compared to E&i, the diffusion barrier then is the minimum value 
of all E&\ of the possible paths that connect R a a and Rad', and the lowest 
energy saddle point is called the transition state. Although often only the 
path with the most favorable energy barrier is important, it may happen 
that several paths exist with comparable barriers or that the PES consists 
of more than one sheet (e.g. Ref. |23j]). Then the effective barrier measured 
in an experiment or a molecular dynamics (MD) simulation represents a 
proper average over all possible pathways. 

In the previous description it was assumed that an adatom moves from 
one binding site to the nearest neighbor one. However, at higher temper- 
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atures diffusing adatoms may from time to time jump over long distances, 
spanning several lattice spacing |Z5[1 . Only little is known about this pro- 



cess. In a recent experimental work on the diffusion of Pd on W (211) Senft 
and Ehrlich j2§j] have extracted from their field ion microscopy (FIM) mea- 



surements an activation barrier for long jumps roughly twice that for single 
jumps. From the analysis of their experimental data they have determined 
the temperature dependent probability for the occurrence of very long jump 
(at least three nearest neighbor distances). These values differ of at least 
one order of magnitude from the theoretical ones [^7j and the reason for 



this discrepancy are still unknown. More effort has to be put into a bet- 
ter understanding of the influence of such long jumps on the intralayer 
transport. 

But diffusion might also occur with a completely different mechanism, 
the so-called diffusion by atomic exchange (or exchange mechanism). The 
adatom can replace a surface atom and the replaced atom then assumes 
an adsorption site. This was first discussed by Bassett and Webber Q8j and 
Wrigley and Ehrlich @. Even for the crystal bulk, namely Si, exchange 
diffusion has been discussed [|28j . This mechanism is activated by the desire 



of the system to keep the number of cut bonds low along the diffusion 
pathway. On fee (100) surfaces diffusion by atomic exchange was observed 
and analyzed for Pt |TT[] and Ir jEfl . For Al(100) it was predicted by 



Feibelman [[T3j and for Au (100) by Yu and Schefher [ptfl . The geometries 
for hopping and exchange diffusion at a fee (100) surface are shown in 
Figs. 0(a) and (b). 

Diffusion along a step edge can also occur via the exchange mechanism as 
illustrated in Fig. 0(c) for a {111} step on the (111) surface. An adatom at 
this step edge experiences a rather high diffusion barrier, if the mechanism 
would be hopping: Either it has to move ontop of a substrate atom or to 
leave the step edge to reach an adjacent step edge position. However, if an 
in-step atom would move out of the step and the adatom fill the opened 
site, the coordination of all the particles would not decrease appreciably 
during the whole process. Thus, the corresponding energy barrier may be 
lower than that of the hopping process. It is also plausible that in general 
the attempt frequency is different for a hopping and an exchange process. 

So far we have only discussed growth processes in the submonolayer 
regime. But except when coverage is sufficiently low atoms might also 
land on top of an existing island. Several questions arise at this point. If 
the island is large enough and the adatom is far enough away from the edge 
of the island, diffusion of this adatom usually will not differ from diffusion 
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Figure 6. Schematic representation of the potential energy close to a step: E d is the 
step-edge barrier, whereas E l d and E% are the diffusion barriers at the upper and lower 
terrace. The additional step-edge barrier is E d 
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on the flat terrace. However, this assumption is not always valid. The 
strain present on the island may affect the self-diffusion barrier. Moreover, 
the atomic structure ontop of the island may differ from the structure of 
the flat surface. An example is given by Pt islands on Pt (111): at T = 
640 K STM images show reconstructed and unreconstructed terraces in 
coexistence with different island densities [|29fl . Furthermore, what happens 
if the adatom is close to the island edge? Is the atom attracted by the edge? 
Does it stay on top of the island or is it hopping down? 

It has been found first by Ehrlich and Hudda |3l| and Schwoebel and 
Shipsey [fnj] and afterwards by a number of other studies that metallic sys- 
tems often exhibit an additional barrier hindering the diffusion over a step 
edge as it is illustrated in Fig. g. This step-edge barrier is often referred to 
as Ehrlich-Schwoebel barrier. Intuitively its occurrence can be understood 
if hopping is the relevant mechanism by employing simple bond counting 
arguments: The atom that diffuses over the step edge is weaker bound right 
over the step edge because at this position the number of bonds is reduced. 
This argument is valid for a step down process by hopping. The situation 
may be different in the case of the exchange mechanism (cf. Fig. 0) because 
here the number of cut bonds remains low along the diffusion pathway. For 
some metallic systems [for example, Al on Al(lll) and Ag on Ag(100)] 
calculations have shown that this is the favored situation (see Sections [§7T] 



13 






Figure 7. The motion of an atom from the upper terrace to the lower terrace down by 
the exchange mechanism. 



and |OQ . Note that the description of the process in terms of a PES [see 
Eq. (|5|)] valid for simple jumps of an adatom (i.e., diffusion by hopping) 
holds for the diffusion by atomic exchange as well [see Fig. §(b) and (c)]. 
The step-edge barrier determines in homoepitaxy whether the growth 
mode is three-dimensional island growth or two-dimensional layer- by- layer 
growth. From experiments this barrier has been estimated by analyzing 
STM images [j32[ , |53fl . The idea is to measure the size of an island just 
when nucleation on top of an island starts and to utilize nucleation theory 
to estimate the step-edge barrier. Similar, indirect studies that interpret 
Monte Carlo simulations [fT3j will give an estimate of the step-edge barrier. 



However, these approaches do not distinguish between different step types 
and are unable to identify the microscopic mechanism for interlayer mass 
transport. 

2.2. Rate equations 

Processes (a), (6), (c), (d), and (e) of Fig. @ form the basis of phenomeno- 
logical rate equations of the form 

^ = $ - 2/^Aq 2 - Aq £ Kj Nj + 2 l2 N 2 + £ 1] N ] (6) 

at j>i j>2 

dN- 

-^ = N 1 ( Kj - 1 N j - 1 - Kj N j )->y j N j + >y j+1 N j+1 . (7) 

These equations describe the time evolution of the adatom density, Aq, 
and the density of islands of size j, Nj : for growth on a flat surface in the 
submonolayer regime. Adatoms are deposited onto the substrate at a rate 
$ = F J\f adatom/s where F is the flux in ML/s and J\f the number of 
atoms pro ML. The second and third term in Eq. (|6]) account for isolated 
adatoms being "lost" because two adatoms can meet at a rate k,\ to form 
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a new nucleus, or adatoms get captured at a rate Kj by an island of size j. 
The last two terms in Eq. (^j) describe further supply sources of adatoms 
and are gain terms because dimers may dissociate and adatoms detach 
from an island of size j at a rate jj. Equation (0) reflects the fact that 
the number of islands of size j increases because islands of size j — 1 grow 
and islands of size J + 1 shrink. The number of islands of size j decreases 
when islands of size j either shrink or grow. Note that no evaporation 
into the gas phase is included in Eqs. (§) and (0) (that means that the 
description in terms of Eqs. (0) and (0) is appropriate only at not too high 
temperatures). One could assume that the rate coefficients Kj oc D and 
7j are independent on the island size and surface coverage (point island 
models) . Another plausible choice is that the rates depend on the length of 
the perimeter of the island so that a first approximation is k,j oc 7j oc y/j 
for compact islands. It has been shown |55j] that the dependence on size 
and coverage is more complex but this will not be discussed here. The rate 
coefficients k s and 7 S are only effective parameters and the physics behind 
them remains unclear. For example, we will see below that the different 
processes shown in Fig. may have different energy barriers and different 
pre-exponential factors, and/or they may proceed by different microscopic 
mechanisms (see Section 0). In principle, all these features can be taken 
into account in Eqs. (0) and (|7j) by introducing different coefficients Kj and 
7j but they become less tractable and their clarity gets lost. A rate equa- 
tion analysis may help to gain some qualitative understanding of growth 
processes but one can not expect insight into the microscopic mechanisms 
governing growth. 

With the assumption that agglomerates of i* + 1 and more adatoms 
are stable against break-up (7^ = for j > i*) one derives the scaling 
relation Qgg| 




i*+2 



iV ls oc — (8) 



where A^ 1S = J2j>i* Nj is the island density and D the diffusion coefficient of 
an adatom on the flat surface. The number i* is called the size of the crit- 
ical nucleus. Relation (@) can be used to extract microscopic parameters 
from experimental (in particular STM) measurements: If one measures the 
island density as a function of F, one can determine the critical nucleus 
size i*. With a known value for i* (assuming that it does not change with 
temperature) one can determine the diffusion barrier E& and the prefac- 
tor Dq if one measures the temperature dependence of A^ 1S . It is not the 
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purpose of this article to review scaling theory; however, we would like to 
point out that this method is not as straightforward as often believed. For 
example, as pointed out by Ratsch et al. [07| , |55f| the size of the critical 
nucleus i* is not always well defined unless the temperature is sufficiently 
low and i* — 1. 

2.3. Critical island area and the action of surfactants 

The definition of layer-by-layer growth mentioned at the beginning of 
Section [I] trivially translates into the following equation [j39|-|r[j] , 



A i* > i/jv is , (9) 

where A 1 ^ is the island area at which nucleation sets in on top of an already 
existing island, and iV ls is the island density [cf. also Eq. (|8j)]. As before 
and in all what follows we assume that we are in a regime of island growth 
rather than step flow. If Eq. (0) is fulfilled, the islands will coalesce before 
a second layer has started to grow. The island density, iV ls , and the critical 
island area, AJ? are controlled by the growth conditions (deposition rate and 
temperature) as well as by the different energy barriers and interactions 
of the deposited adatoms on the surface and by the minimal size of an 
island nucleus. A^ is determined by the probability that a number of i* + 1 
adatoms meet on the same island and form a stable nucleus. 

Without a step-edge barrier the adatoms that land on an island are not 
hindered to move down and bind at the favorable sites. Thus, the formation 
of a nucleus on top of an island becomes unlikely, and layer-by-layer growth 
is expected. However, the situation is different when atoms that land on 
top of an island are hindered by a step-edge barrier to move downwards. 
In this case it is more likely that i* + 1 adatoms meet to form a stable 
cluster that subsequently will grow into a bigger and bigger island. Thus, 
A 1 * is smaller (compared to l/iV ls ) and it is more likely that islands on the 
surface reach the area A 1 * before the layer is completed. In other words, 
when a noticeable step-edge barrier exists Eq. (0) may not be fulfilled and 
the system will grow three-dimensionally 

Eq. (|9D can be also rewritten in the following form 



6 C = AjfiV 18 > 1 , (10) 

where G c is the critical coverage for which, when it is exceeded, islands will 
grow on top of already existing islands. Looking at experimental situations 
it appears that the conditions set by Eqs. (|9j) and (pD|) are slightly too 
strong and it is probably sufficient to request O c > 0.9 for good layer-by- 
layer growth. The importance of the above equations is that they show 
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that the growth mode can be influenced in two independent ways: One can 
modify the island density (which is controlled by the adsorbate mobility 
at flat regions) or one can modify the critical island area A™ which is 
largely determined by the physics at step edges. Thus, it follows that the 
growth mode can be changed from three-dimensional to layer-by-layer in 
the following ways: 

1) It has been shown by Kunkel et al. [f|2j] that on Pt (111) at very low 



temperatures the island shape is fractal. This implies that the island 
perimeter is particularly long, that the islands have rough edges, and 
many rather thin branches. As a consequence the step-edge barrier 
might be reduced substantially. But even if the step-edge barrier 
remained unchanged, the probability for an adatom which lands on 
such an island to move down to the lower terrace is high, because 
the adatom will visit the edge very frequently. This makes it likely 
that a possibly existing energy barrier can be overcome. Moreover, a 
large number of kink sites or weakly bound atoms is present along the 
edges of the fractally shaped islands. Thus, an exchange downward 
diffusion of an adatom on the upper terrace at these sites may be 
more likely than the same mechanism involving an atom of a compact 
step edge. Furthermore, we note that the island density iV ls is high at 
low temperature [cf. Eq. (§)] which reduces the probability that more 
than one atom land in a reasonable time interval on the same island. 
Altogether these properties give rise to layer-by-layer growth. 



2) A second possibility was demonstrated by Rosenfeld et al. Qffil . These 



authors have shown that increasing the island density A^ 1S is in fact 
sufficient to achieve layer-by-layer growth. A high island density can 
be obtained for example by lowering the temperature or increasing 
the deposition rate in the very beginning of growth. Once the island 
density is increased, the growth is continued at normal (higher T and 
lower F) conditions. Thus, the island density A^ 1S is set large (to low 
T and high F parameters) but A™ is not reduced and remains at the 
value determined by the "normal" T and F parameters. Indeed, three- 
dimensional growth did not start before the layer was completed. 

3) For completeness we mention the possibility to enhance the mobility 
of deposited adatoms by photo-stimulation. However, we will not 
elaborate on this mechanism. 

4) A very interesting way to achieve layer-by-layer growth uses surface 
contaminants, so called surfactants. There is one necessary condition 
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these species should fulfill: They should stay on the surface during 
growth, thus they should not become buried during the growth pro- 
cess. While a good probability of surface segregation is necessary it 
alone would not affect the growth mode. There are the following pos- 
sible mechanisms that, when active, provide that a surface segregating 
contaminant increases the inter-layer mass transport: 

i) The simplest idea that comes to mind is that the surfactant deco- 
rates edges of steps and islands and reduces the step-edge barrier, 
since the atom- surfactant interaction is usually weaker than the 
atom-atom interaction. Figure demonstrates how this could be 
achieved. Lowering the step edge barrier facilitates the interlayer 
transport and O c « 1. Oxygen for Pt(lll) j33| and indium for 
Cu(100) seem to have this effect |55j and enhance the 2D char- 



acter of the growth mode. Recently another picture has been 
proposed |E6|] : A repelling action of In at steps hampers the at- 
tachment of Cu atoms approaching the step edge from the lower 
terrace and gives rise to an enhancement of the island density. 
However, the experimental data of Ref. |35j seems to rule out this 
scenario. 

ii) It is also possible that surface impurities induce a potential en- 
ergy gradient that attracts deposited atoms towards the step; for 
deposited atoms that land on an island the number of visits at 
the edge is thus increased, and as a consequence the probability 
to move down is increased as well. 

Hi) Surfactants may act as nucleation centers, thus increasing the 
island density iV ls , while A 1 * remains unchanged (of course, an 
increase would be even better). This will induce layer-by-layer 
growth, provided that the probability that atoms which land on 
an island and move to the lower terrace is not reduced as well. 
Moreover, if the surfactant increases the diffusion barriers E l d and 
E% but leaves E^ (cf. Fig. |6j) essentially unaffected, it follows 
that E^ — E% is reduced and the wanted effect may result. This 
mechanism was recently discussed by Zhang and Lagally [||7| . 



iv) A forth possibility was discussed in the context of the surfactant 
action of Sb on Ag(lll). The basic mechanism here is that Sb 
impurity atoms on the surface are practically immobile and act 
repulsively to deposited Ag adatoms. This will also increase the 
island density A^ 1S and thus further two-dimensional growth. This 
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mechanism will be discussed in Section |4.2.2| . 

Tersoff et al. ||l| have recently discussed Eq. (0) by assuming a circu- 
lar island shape and various sizes for the island nucleus. They demon- 
strated that the critical island area is indeed a useful concept. However, 
we hesitate to give an explicit formula for it in terms of diffusivities and 
deposition rate because in reality A™ will depend sensitively on the long- 
and medium-range adatom-adatom and adatom-step interactions (see the 
attractive gradient towards the step on the upper and lower terraces in 
Fig. P), as well as on the diffusion barriers of adatoms parallel to step 
edges, as these determine the actual shape of an island. 

A study of Memmel and Bertel [£|8| has raised an interesting point. They 



propose a simple model which connects the diffusion behavior on metal 
surfaces to the charge density supplied by occupied two-dimensional free- 
electron surface states. The argument is very appealing: A decrease in 
the difference between the step edge barrier and the activation energy for 
diffusion on the flat terraces could enhance the interlayer mass transport. 
The barrier for the diffusion on the flat terrace is mainly determined by 
the corrugation of the electron density to which both bulk Bloch states 
as well as surface states contribute. The surface states are particularly 
interesting, since they can strongly be influenced. A depopulation of these 
states induced by confinement onto small islands or by the presence of an 
appropriate surfactant increases the diffusion barrier on the flat surface 
with a consequent reduction of the additional step edge barrier at the step 
edge. Thus, an increased interlayer transport is expected with the related 
layer-by-layer growth. This picture seems to be appropriate for the effects 
of oxygen on the growth mode of Pt on Pt (111) [f|3| . 



3. TOTAL ENERGY AND THE DESCRIPTION OF GROWTH 

In Section gTT] we defined the potential energy surface (PES) of a diffus- 
ing adatom. Obviously, the PES is governed by the breaking and making 
of chemical bonds, and we also noted the need to take atomic relaxations 
into account [cf. Eq. (0)]. Thus, the evaluation of the PES requires an 
accurate, quantum-mechanical description of the many-electron system. 
This can be achieved by modern density functional theory calculations 
that combine electronic self-consistency and efficient geometry optimiza- 
tion. Approximate methods, based on the concepts of DFT but employing 
empirical parameter instead of elaborate calculations have been developed 
as well. Such approximate methods are widely used by several groups to 
investigate surface properties and to perform MD investigations of adatom 
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diffusion. We will sketch their main characteristics in Section PTT] . A de- 
scription of the basic concepts of DFT is then given in Section |5T2| , and 
Section |3.3| describes how DFT is implemented into accurate self-consistent 



calculation methods. Finally, in Section |3.4j we describe briefly the kinetic 
Monte Carlo (KMC) technique that is capable to tackle the realistic time 
and length scale of growth. 

3.1. Bond-cutting methods 

Several methods have been developed based on the idea that the en- 
ergy of a many-electron, poly-atomic system can be written in terms of 
contributions from the individual atoms: 

E tot ({R I }) = Y,E I . (11) 

The sum goes over all atoms, and Ej is the contribution of the 7-th 
atom. Ei depends sensitively on the local geometry of atom I (its em- 
bedding). The different bond-cutting methods differ in the treatment of 
the actual form of the "embedding function" and in the way to determine 
the necessary materials parameters. The differences are not very signifi- 
cant and the most popular names of these methods are: embedded atom 
method (EAM) ggJEgJ , effective medium theory (EMT) |5IH53fl, Finnis- 
Sinclair A-body potentials j55fl , second-moment approximation [|55fl , and 
glue-model Q57[ 



In the simplest version of a bond-cutting approach it is assumed that the 
energy per atom Ei varies linearly with the atom's coordination number. 
Thus, it is assumed that the strength of a bond is invariant of the number 
of bonds the atom does form. This approach clearly neglects the quantum- 
mechanical properties of bonding, namely that the bond strength saturates 
at a certain number of neighbors [|5B| , |5B| , |5^|1 . In fact, detailed DFT studies 
have shown [|5B],|59j that the dependence of Ej on the local coordination is 



very similar to 

E^-A^fCi + BCi , (12) 

with Cj the coordination number of the 7-th atom. 
A more general approach gives 

£/ = F(p/) + i£0(|R.j-R.i|) • (13) 

Here (j) describes the pair-wise, repulsive interaction between atoms, and 
F is called the embedding function, that depends on the electron density 
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created at site /. For the effective medium theory Christensen and Ja- 
cobsen |5Dj have shown that Eq. (P~3| ) resembles the behavior of the simple 
function noted in Eq. (|T21), but also contains some refinements. Indeed, is 
has been shown that Eqs. (|T1| ) and ([13]) represent an approximation of the 



total-energy expression of density functional theory []53|,|53j] 

The main problem in actual calculations is to determine the necessary 
parameters to define the embedding function. Typically the parameters 
are obtained by fitting results from a treatment based on Eq. ( P~5| ) to some 
experimental or DFT results of "related systems". The results depend 
on what systems and what properties are chosen. The predictive power 
of these methods has to be questioned (see e.g. [pl| , |62|1 ) . We note, how- 



ever, that bond-cutting methods hold a significant share of the quantum- 
mechanical description and thus are most valuable to summarize and to 
explain trends of results obtained by DFT calculations. 

3.2. Density functional theory 

The total energy of an iV e -electron, poly-atomic system is given by the 
expectation value of the many-particle Hamiltonian, using the many-body 
wave-function of the electronic ground state. For a solid or a surface the cal- 
culation of such expectation value is impossible when using a wave- function 
approach. However, as has been shown by Hohenberg and Kohn [ |53f| , the 
ground-state total energy can also be obtained without explicit knowledge 
of the many-electron wave-function, but from minimizing an energy func- 
tional E[ri\. This is the essence of density functional theory (DFT), which 
is primarily (though in principle not exclusively) a theory of the electronic 
ground state, couched in terms of the electron density n(r) instead of the 
many-electron wave function \I/({rj}) with r$ the coordinates of the i-ih 
electron. 

The important theorem of Hohenberg and Kohn |53j (see also Levy |j53f ) 



tells: The specification of a ground state density n(r) determines the many- 
body wave function. In other words, Hohenberg and Kohn realized that for 
the ground state the known functional n(r) = n[\P] = (\I/| J2iS(r — r^) |M/) 
can be inverted, i.e., \I/ = \£[n(r)]. Although it was shown that \£[n] exists, 
its explicit form remains unknown. 

With the help of this theorem the variational problem of the many- 
particle Schrodinger equation transforms into a variational problem of an 
energy functional: 

E < {^\H C \^) = E v [^[n}} = E v [n] . (14) 

Here Eq is the energy of the ground state, H e the Hamiltonian of the elec- 
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trons, E v [n] = / dr v ext (r)n(r) + G[n], and v ext (r) is the external potential 
(typically ?; ext (r) is the Coulomb potential due to nuclei). In this functional 
n(r) is the variable (the electron ground-state density of any iV e -electron 
system), and v ext (r) is kept fixed. G[n] is a universal functional indepen- 
dent of the system, i.e., independent of v cxt (r). For example, G[n] is the 
same for an H-atom, a CO-molecule, a solid etc.. The main advantage of 
this approach is that n(r) only depends on three variables while ^ ({r^}) de- 
pends on many variables (the 3iV e coordinates of all electrons) J65|. Thus, 
it is plausible that the variational problem of E v [n] is easier to solve than 
that of (^\H e \ty), yet the result for the ground-state energy and the ground 
state electron density will be the same. The total energy entering Eq 
is 



E M ({Rj}) = E ({Rj}) + \ £ ZjZ * , (15) 

where {Rj} includes all atoms, and Zj is the nuclear charge. 

Earlier work (in particular the Thomas-Fermi approach) had shown that 
the treatment of the kinetic energy (^| — |V 2 |^) is of particular importance 
and Kohn and Sham |j67[] therefore wrote the energy functional in the form 



E v [n] = T s [n] + J dr v cxt (r)n(r) + - J dr v R (r)n(r) + E xc [n] , (16) 

where T s [n] is the functional of the kinetic energy of a system of non- 
interacting electrons with density n(r), and -u H (r) = J dr' rfj, , the Hartree 
potential, is the time-averaged electrostatic potential created by the elec- 
tron density, and E xc [n] is the so-called exchange-correlation functional. It 
accounts for the Pauli principle, dynamical correlations due to the Coulomb 
repulsion, and the correction of the self-interaction included for convenience 
in the Hartree term. 

With Eq. ( \T6j ) the problem of the unknown functional G[n] is transformed 
to one involving T s [n] and E xc [n] . The functional defined by Eq. (|T5|) can be 
also modified by adding terms which vanish at the correct electron density. 
Such a new functional E v [n] may converge faster towards the ground state 
or may depend less sensitive on the input density. The latter implies that 
the input density does not need to be very accurate, yet the resulting 
energy represents an acceptable approximation for the correct total energy 
(see e.g. Ref. [|68f| ). Although the functional T s [n] is not known explicitly 
in a mathematically closed form, it can be evaluated exactly by using the 
following "detour" proposed by Kohn and Sham. The variational principle 
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applied to Eq. ( |nj|) leads to 



Sn(r) 5n(r) l ' Sn(r) 

+ v cS (r)=/i , (18) 



MM , .c ff , 



Sn(r) 

where \i is the Lagrange multiplier associated with the requirement of a 
constant particle number and thus equals the electron chemical potential. 
The effective potential is defined as 

v cS (r) = v c *\r) + v R {r) + v* c {r) , (19) 

with v xc (r) = 5E xc [n]/5n(r), and n(r) is a ground-state density of any 
non- interacting electron system, i.e., 

n(r) = £/«l&Mr , (20) 

where we introduced the occupation numbers fa which is simply the Fermi 
function. Because T s [n] is the kinetic-energy functional of non-interacting 
electrons, Eq. fllgp [together with Eq. (|2D|) 1 is solved by 



--V 2 + v cS ( 



(21) 



These are the Kohn-Sham equations. They have to be solved self-consistently 
together with Eqs. ((n|) and (|2D|). In principle this gives the exact ground- 
state electron density and total energy of a system of interacting electrons. 
The above Eqs. ( |I5| ) - (|2T| ) contain one term which is not known exactly. 
This is the exchange-correlation functional £" xc . For a better understand- 
ing of this term it is instructive to introduce the following Schrodinger 
equation [|5gH7If| : 

Hx\*x) = {T + V x + \W}\* x ) , (22) 

where A defines the strength of the electron-electron interaction: A = 
corresponds to the noninteracting system and A = 1 to the interacting 
physical system. V\ — J2i v\(ri) is an external potential chosen to maintain 
the ground state density n(r) at its A = 1 value independently of A. The 
exchange-correlation energy E xc can be written in a form which resembles 
a Coulomb interaction, 

E xc [n] = I I [ drdr' n(r) , l , n xc (r, r') , (23) 

2 J J r — r 
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Figure 8. Sketch of the pair correlation functions g^ and g^ of a homogeneous electron 
gas in the Hartree (dotted line) and Hartree-Fock (dashed line) approximation, and in 
the exact treatment (DFT) (solid line). 



where n xc (r, r') is called the exchange-correlation hole. It can be written 



as 



n xc (r,r'1 



n(r')[£(ry)-l] 



(24) 



with the electron pair-correlation function g(r, r'). Using a wave function 
formalism one then obtains []7Tj 



g(r, r') = j Q dX 



[Vx\n{r)n{r')\yx) S(t - r 1 ] 



n[T)n[r 



n[T] 



(25) 



where n is the density operator. According to Eq. ( |23|) , the exchange- 
correlation energy arises from the Coulomb interaction of each electron 
(e.g. the one at r) with a charge distribution n xc (r, r'), i.e., the exchange- 
correlation hole surrounding that electron. The hole is a consequence of 
the exchange and Coulomb interactions that cause a depletion of electron 
density in the vicinity of each electron. We note that the picture behind 
the depletion is a dynamical one. In the time average the depletion is not 
seen in the electron density, but it gives rise to an important lowering of 
the total energy, compared to that of a non-interacting electron system. 

It is instructive to compare Hartree, Hartree-Fock, and DFT descriptions 
and in Fig. |8] we show schematically the spin-resolved pair-correlation func- 
tions g^i and g^ of a many-electron system of constant density (also called 
a jellium system). The quantities g^ and g^ are the pair-correlation func- 
tions for electrons with parallel and antiparallel spins, respectively. In a 
spin unrestricted calculation g of Eq. (|Z5|) is given by g = g^ + g^. In 
the pure Hartree description exchange and correlation are ignored and the 
pair-correlation functions are constant and equal to 1/2. The Hartree-Fock 
method accounts for dynamical correlations due to the Pauli principle (the 
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exchange). The pair-correlation function g^ : which is the probability of 
finding an electron at r' and an electron at r with the same spin, shows a 
dependence on the interelectronic spacing: The closer we are to the elec- 
tron at r, the lower is the probability to find another electron. However, g^ 
is constant because the Hartree-Fock description only contains the "Pauli 
correlation" which affects electrons with the same spin, but the depletion 
due to Coulomb repulsion which is independent of the spin is neglected. 
These are accounted for in DFT. Thus, DFT gives a correct description of 
the fact that electrons move in a correlated way and that this correlation 
is due to the Pauli repulsion (for electrons of equal spin) and the Coulomb 
repulsion (for all electrons). 

The problem that remains in an actual DFT calculation is that the func- 
tional E xc [n] is unknown. However, some general properties of this func- 
tional and values for some special cases are known. Detailed and accurate 
understanding exists for systems of constant electron density. The asymp- 
totic behavior at low and high densities is given by expressions derived by 
Wigner |72]] and Gell-Mann and Brueckner |73j and for intermediate den- 



sities quantum Monte Carlo calculations have been performed by Ceperley 
and Alder |73j] . This gives the simple curve shown in Fig. 0, and this result 



for e xc (r s ) = Clda! 71 ) * s then used in the functional 

El^in] = J dr n(r)e^ A (n(r)) , (26) 

which is the local-density approximation (LDA) |T7|. Thus, in the LDA 
the many-body effects are included such that for a homogeneous electron 
gas the treatment is exact and for an inhomogeneous system exchange 
and correlation are treated by assuming that the system is composed from 
many small systems with a locally constant density. 

The LDA can be improved by including the dependence on the density 
gradient which leads to the generalized gradient approximation (GGA). 
Several different GGAs were proposed in the literature [|75H50fl and have 
been used successfully for DFT calculations of atoms, molecules, bulk 
solids, and surfaces (an overview can be found in Refs. [jSOUSIfl ), but also 



limitations have been identified for example by Mitas et al. j82f and Um 



rigar and coworkers [p3fl . It is by now clear that the lattice constants 



calculated with a GGA are typically larger than those obtained with the 
LDA, with the experimental values usually being in between. Binding 
energies (or cohesive energies) of molecules and solids as well as energy 
barriers of chemical reactions are improved by the GGA (see Ref. [|S^ , |S5 



and references therein). Still, for surface diffusion DFT-LDA calculations 
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Figure 9. Exchange-correlation energy per particle, e xc , of homogeneous electron gases 
with density parameters r s . The electron density and the density parameter are related 
by n 
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give energy barriers in good agreement with those deduced from experi- 
ments and with GGA calculations. Thus, although the total energies are 
changed when going from the LDA to the GGA, the changes in energy 
barriers that are the differences between total energies are typically less 
pronounced (see e.g. Ref. |j62|J85fl ). 



3.3. Implementation of DFT into state-of-the-art computations 

Typically there are only a few possible adsorption sites and channels 
for diffusion. This is illustrated in Fig. [TU| for the fee (111) and fee (100) 
surfaces. For adatoms that are chemically similar (or equivalent) to those 
of the substrate the stable sites are those with high coordination and for 
hopping diffusion the transition state is at the bridge site. The relevant 
information about the PES then is obtained by calculating the total energy 
of the system with the adatom placed in those positions. In general, more 
care is necessary because the bridge site could also be a local minimum of 
the PES and the energy barrier could be in between the high coordination 
and the bridge sites. Furthermore, it is possible that the diffusion does not 
proceed by hopping but by atomic exchange [|SKHfl. 

In the bulk crystal the three-dimensional periodicity can be exploited 
by using Bloch's theorem. Unfortunately, the presence of a surface and 
an adatom on top of it breaks all periodicities. The (in principle) best 
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Figure 10. Top view at a fee (111) and (100) surface. The adsorption sites labeled fee, 
hep, and hollow site usually correspond to the most stable binding sites while the bridge 
site is the transition state of a hopping diffusion. 



approach to treat such difficult situation is given by the Green-function 
method |S7| , |S5fl . An, at least in the past, popular approximation for ad- 



sorbate systems is the cluster approach |]89||. The presently most efficient 



and practical approach that was also used in the results discussed below is 
the supercell approach. The supercell may also be called a big cluster, but 
in contrast to conventional cluster calculations the supercell is periodically 
repeated. As a consequence, the cluster boundary is treated physically 
very accurately, and by utilizing the periodicity it is possible to use very 
big cells. The idea of an adatom on top of a substrate in the supercell 
approach is sketched in Fig. [IT]. The adatom is placed on top of a slab 
of a certain number of layers. The number of layers Af s i a b has to be suf- 
ficiently large so that the adatom does not "feel" the presence of another 
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Sketch of a supercell describing an "isolated" adatom at a surface (side 



surface on the other side of the slab (or at least that the quantity to be 
computed, such as a diffusion barrier, is not affected by the other surface). 
Alternatively, one could also place an adatom on either side of the slab; 
in this case, there are more symmetries in the geometry but more layers 
are needed in the slab to screen the mutual interaction between the two 
adatoms through the slab. The adequate number of layers iV s i a b depends 
on the properties that one wants to calculate and the surface orientation, 
and careful tests have to be carried out. For example for the Ag (111) sur- 
face four layers are sufficient when the adatom is placed on only one side 
of the slab, while for Al (100) seven layers are necessary. 

As illustrated in Fig. |TT] the geometry repeats periodically in vertical 
and lateral directions. The lateral periodicity implies that a single adatom 
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placed on a substrate is not at all a single adatom; if the cell size parallel 
to the surface is chosen for example as (2 x 2) we actually calculate a 
system with a coverage of 25 %. It is therefore important to test that 
the interaction with the neighboring adatoms can be neglected. On a 
(111) surface a cell size of (2 x 2) is usually sufficient, but sometimes 
larger cells [(3 x 3) or even (4 x 4)] are necessary. To model a diffusion 
event along or across a step one either chooses a small island on top of 
a substrate or a vicinal surface. The latter has the advantage that there 
is only one step in the unit cell so that a smaller cell size is required in 
order to attain a negligible step-step interaction. The system also repeats 
in vertical direction separated by a vacuum region. The thickness of the 
vacuum region has to be tested as well, but the computational cost of a 
thicker vacuum region is relatively small compared to a larger cell size 
or a higher number of slab layers (for a deeper discussion of the above 
technicalities see e.g. Ref. 0). 

Core electrons typically do not take part directly in the binding process 
of atoms in molecules and solids, and the nature of the chemical bond is 
mainly determined by the valence electrons. This is exploited by the frozen 
core approximation where the core electrons are effectively combined with 
the nuclei to form frozen unpolarizable ions. Still, not just the electrostatic 
potential but also the quantum nature of the core electrons is felt by the 
valence electrons. For example, different wave functions have to be orthog- 
onal and therefore the valence wave functions have nodes and oscillate in 
the core region. For practical calculations one needs to expand the wave 
function in a suitable basis and we choose a plane wave basis set |JT] 



fe(k,r)) = £c„ k (G)|G + k) . (27) 

G 

A plane wave description of wave functions that have nodes and oscillate 
requires a very large number of plane waves. This inconvenience is cured 
efficiently by the pseudopotential approach. Modern ab initio pseudopo- 
tentials reproduce the potential of an atom exactly outside the core region 
defined by a radius r c and are rather smooth inside the core region. An im- 
portant requirement on a "good" pseudopotential is that it is transferable. 
This means that the pseudopotential should behave like the all-electron 
potential in a variety of different chemical situations. Pseudopotentials 
that reproduce the same charge inside the core region as the all-electron 
potential, and therefore have the same scattering properties, are referred 
to as norm- conserving. Those that are often used have been developed 
by Bachelet, Hamann, and Schliiter |J2]], Troullier and Martins [|93f| , and 
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Gonze, Stumpf, and Scheffler |j9^fl . Recently, Vanderbilt [£55| proposed ab 
initio pseudopotentials that drop the condition of norm-conservation and 
therefore can be used with a lower number of plane waves. The gain in com- 
puter time due to the smaller basis set is partially compensated by the costs 
to calculate the correction required by the neglected norm-conservation. 
The electron density is calculated according to Eq. (^0|) as 

«(r) = EE^(k))|«k,r)| 2 (28) 

k j 

where the integration over the Brillouin zone has been replaced by a sum 
over a mesh of k-points and Wk is the weight of the k-point. A conve- 
nient scheme to construct an appropriate k-point mesh is described by 
Monkhorst and Pack [|96|| . In ab initio pseudopotential calculations some 



matrix elements and some integrals are efficiently evaluated in real space 
whereas others are efficiently evaluated in reciprocal space. The technique 
of fast Fourier transformation enables a numerically fast change from one 
representation to the other. Technical details of the computational proce- 
dures are described for example in [|97H99(| • 

3.4. Kinetic Monte Carlo approach 

MD simulations can provide important insight into elementary micro- 
scopic mechanisms but typically they cannot be used for growth studies. 
The time between two successful diffusion events is often of the order of 
nanoseconds. During this time the adatoms undergo several (e.g. 10 3 and 
more) unsuccessful attempts. Since MD calculates all these unsuccess- 
ful atomic movements explicitly, MD simulations can cover times of some 
picoseconds, possibly some nanoseconds and therefore are usually inappro- 
priate to describe the spatial and temporal evolution of growth patterns, 
that typically develop on a time scale of seconds. Instead, the method of 
choice for studying the spatial and temporal development of growth is ki- 
netic Monte Carlo (KMC) gUtJ. The key idea behind KMC is to describe 
stochastic processes (such as deposition, diffusion, desorption, etc.) on the 

microscopic scales by rates T& = Iq exp(— E^' /k^T), that were discussed 
in Section [271] above. In KMC a process j is thermally activated with a 
relative probability given by: 

w> = — = ^eM-E% ) /k B T) . (29) 

We define the total rate R = £j r^, where the sum runs over all possible 
processes. Deposition is accounted for in this description by the deposition 
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rate Tq — F. Equation (|29|) satisfies the condition of detailed balance. Us- 
ing a stochastic approach the explicit calculation of unsuccessful attempts 
is avoided. Yet, the result of a KMC study will be the same as that of a MD 
simulation, provided that the underlying PES is the same. The strategy 
of KMC can be summarized as follows: 

1) Determination of all processes j that possibly could take place with 
the actual configuration of the system. 

2) Calculation of R. 

3) Choose two random numbers p\ and pi in the range (0, 1]. 

4) Find the integer number / for which 

£ r w <piR<J2 r0) • (so) 

3=0 j=0 

5) Execute process /. 

6) Update the simulation time t := t + At with At = —In(p2)/R . 

7) Go back to step 1). 

Step 6) ensures that a direct and unambiguous relationship between KMC 
time and real time is established, since the KMC algorithm effectively 
simulates Poisson processes. 
KMC differs from the algorithm proposed by Metropolis et al. [jlUl 



that was successfully employed to determine the static properties of many- 
particle systems. In the Metropolis Monte Carlo scheme the probability 
that a new configuration is accepted is proportional to exp(— AE/k^T) 
where AE is the difference between the total energies of the system in the 
new and old configuration. Such an algorithm searches for the configura- 
tion corresponding to the minimum of the total energy, and the sequence 
of generated configurations does not correspond to the real time evolution 
of the system. 

KMC simulations have been used to study crystal growth of semiconduc- 
tors (e.g. |TU2j-|TU5j ) and metals (e.g. [jIQ5HI08|| ) . However, most of these 
studies have been based on restrictive approximations. For example, the 
input parameters have been treated as effective parameters determined 
rather indirectly by fitting to experimental quantities, like intensity oscil- 
lations in helium atom scattering (HAS) measurements or reflection high 
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energy electron diffraction (RHEED), or they were obtained from STM 
studies of island densities. Thus, the connection between these parameters 
and the microscopic nature of the processes may be somewhat uncertain. 
Often the correct surface structure was neglected and the simulation was 
done on a simple cubic lattice while the system of interest had an fee or bec 
structure. Despite these approximations such studies have provided quali- 
tative and in some cases also quantitative insight into growth phenomena. 
It is desirable to carry out KMC simulations with the proper geometry and 
microscopically well founded parameters. This has been done for example 
in Refs. [|1U9H111 |1 where semi-empirical calculations such as the embedded 
atom method or effective medium theory (cf. Section p.l|) have been em- 
ployed to evaluate the PES. However, the most accurate, but also most 
elaborate approach to obtain the PES employs DFT as described in Sec- 
tions ^3 and P73| . Results of a KMC study performed with input from such 



DFT calculations are presented in Section |4."1.2| . 

4. RESULTS FOR FCC (111) AND FCC (100) SURFACES 

We now discuss some recent DFT studies of various growth phenomena 
at (111) and (100) surfaces of fee metals, in particular aluminum (Sec- 
tions |0] and [O]) and silver (Sections 4.2 and ffTjj) . Section |4.1.1| presents 



DFT calculations for the self-diffusion at Al(lll), and these results form 
the basis of the ab initio KMC simulations described in Section |4.1.4 Sec- 



tion |4.5| summarizes the understanding of self-diffusion at Al (100), and in 



Section fO| we identify the microscopic processes on the Ag (100) surface. 



Surface stress can affect diffusion and crystal growth quite substantially. 
Experimental investigations have been carried out for the self-diffusion at 
Ag(lll) and at thin silver films on Pt(lll) [|TT2j] . In Section |4T2T1| we 
present results of recent DFT calculations for these systems. Section |4.4.2| 
discusses the influence of surface stress on the diffusion mechanism, showing 
that stress can even change the mechanisms. 

One goal of epitaxial growth is to achieve atomically flat and defect free 
surfaces of specified crystallographic orientation under convenient growth 
conditions, i.e., not too low deposition rates and not too high tempera- 
tures. In many cases this goal can be reached only with the use of surface 
contaminants that act as surfactants [J113| , |114J1 • Section |4.2.2| presents DFT 



results concerned with the surfactant functioning of Sb for the growth of 
Ag(lll). 

4.1. Growth at Al(lll) 
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4.1.1. Microscopic processes 

On the (111) surface of a fee crystal there are two kinds of close packed 
steps, that were already discussed (see Fig. f|). Although the packing of 
the step-edge atoms is identical in these two steps, the coordination of the 
substrate atoms at the step edge is different: At the {lll}-faceted step 
the substrate atoms have a coordination of 11, and at the {100}-faceted 
step the substrate atom coordination is 10. Thus, a small difference in 
the step-formation energy is expected. If lattice relaxations are negligible, 
the {lll}-faceted step might be slightly favored. Indeed, DFT calcula- 
tions [PD| , |115|1 confirm the trend suggested by the coordination numbers: 
The {111} facet is slightly favored with a^ 111 ^ = 0.232 eV per atom over 
the {100} with a^ 100 ^ = 0.248 eV per atom. Neglecting contributions from 
more open steps, as these have a higher formation energy, one can now 
predict the equilibrium shape of islands. According to the Wulff construc- 
tion in thermodynamic equilibrium the island borders have distances from 
the island center that are proportional to their formation energies. As 
a consequence the equilibrium shape of islands for A1/A1(111) should be 
hexagonal with the edges alternating between those with a shorter {100} 
and a longer {111} facet. The resulting ratio of the lengths L^ 100 ^ / ZA 111 } 
is 4 : 5. 

The activation barriers for the most important processes of self-diffusion 
at Al (100) are collected in Table [T]. In addition we note that the calcula- 
tions predict that an isolated adatom on the flat surface favors the hep site 
(see Fig. |TD|) and that the energies of the bridge and fee positions are almost 
degenerate. Therefore, the diffusion path for an isolated adatom goes from 
one hep site to the adjacent one through fee and bridge positions. Since 
the fee site continues the ABCABC stacking of the fee crystal, whereas the 



Table 1 

Energy barriers E<i for different self- diffusion processes on Al(lll). 

process mechanism E& (eV) 

flat Al(lll) hopping 0.04 

{111} step | exchange 0.42 

{111} step | hopping 0.48 

{100} step | hopping 0.32 

{100} step | exchange 0.44 

{111} step _L descent exchange 0.06 

{111} step _L descent hopping 0.33 

{100} step _L descent exchange 0.08 

{100} step _L descent hopping 0.45 
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equally threefold coordinated hep position belongs to an ABCAC stack- 
ing, this result is somewhat unexpected. The epitaxial continuation of the 
crystal with adatoms occupying fee sites is recovered when the coverage is 
increased (6 > 1/4 ML). 

The DFT calculations also predict a long-range attraction of adatoms 
towards step edges for approach from the upper as well as from the lower 
terrace. It appears that this attraction is actuated by electronic surface 
states. The attraction is weak at long distances but close to the step it 
becomes so strong that particularly at the lower terrace an adatom will be 
funneled toward the step. This is clearly visible in Fig. [T2] (lower curve) 
where the total energy along the adatom diffusion path involving the mi- 
gration toward and over a step and on the flat surface is displayed for the 
hopping (upper curve) and the exchange (lower curve) mechanism. The 
attraction is present with and without relaxation of the system and thus 
cannot be elastic. An electrostatic origin can also be discarded, since the 
dipoles located at the step and of the isolated adatom have the same sign. 
Thus, the resulting dipole-dipole interaction is repulsive. However, the 
adatom and step induce localized electronic states that interact and it was 




Figure 12. Total energy along the diffusion path of an Al adatom over a {lll}-faceted 
step on Al(lll). The upper curve is evaluated for a hopping process, while the lower one 
refers to an exchange process. The generalized coordinate is Q = X\ + X 2 with Xj = 
x-coordinate of the adatom /. The x axis is parallel to the surface and perpendicular to 
the step orientation. For the undistorted step X 2 = 0. 
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concluded [|9(jyil5|1 that they are responsible for the long-range attractive 
interaction. 
Figure [12] also shows the difference between the hopping and the exchange 



mechanism for the diffusion across the step from the higher to the lower 
terrace. The upper curve is calculated for the hopping process, and the 
presence of an energy barrier that hinders the roll over of the adatom from 
the upper to the lower terrace is clearly seen. On the other hand, practically 
no hindrance exists for diffusion by exchange (cf. Fig. 0). The activation 
barrier for the exchange process for the step-down motion is very low and 
comparable to the diffusion barrier on the flat surface: 0.08 eV and 0.06 
eV for {100}- and {lll}-faceted steps, respectively, whereas it is 0.04 eV 
for diffusion on the flat surface. Thus, we predict layer-by-layer growth for 
Al on Al(lll) for a wide range of substrate temperatures. The barriers 
for the exchange might hamper 2D growth only at T ~ 25 K. However, at 
such a low temperature the island edges are frayed which may reduce the 
barriers resulting in layer- by-layer growth. 

One origin for the preference of the exchange process at step edges might 
be the bonding character of Al. Although Al is often considered as a 
jellium-like metal, it is a rather covalent atom with its sp valence electrons 
(we remind that Al and As form the AlAs compound, a covalent, zincblende 
semiconductor). Thus, similarly as discussed by Pandey |28j] for the ex- 



change diffusion in Si bulk, and by Feibelman |T3j] for exchange diffusion at 
Al (100), we believe that this mechanism is favored by the tendency of the 
system to keep a low number of cut bonds along the diffusion pathway. We 
expect that exchange diffusion is a rather common mechanism for down 
movement of adatoms at step edges. 

We now address the diffusion of adatoms parallel to the two close-packed 
steps. As will be shown in Section |4.1.^ , this is of particular importance for 
the shapes of islands which develop during growth. The migration along the 
steps may take place via a hopping or an exchange mechanism. The calcu- 
lations predict that along the {100}-faceted step an adatom preferentially 
jumps to an adjacent site with an activation barrier of 0.32 eV, whereas 
along the {111} facet diffusion by atomic exchange is preferred with ac- 
tivation barrier of 0.42 eV. To understand this difference we consider the 
positions labeled TSjioo} an d TS{im in Fig. [T3| that are the lowest-energy 
transition state of a hop along both steps. Along the {lll}-faceted step 
the adatom in the transition state has three neighbors, whereas in the 
same position along the {100}-faceted step the adatom has four neigh- 
bors |}9U[]. The higher coordination suggests a lower barrier for hopping 
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Figure 13. The transition sites for the hopping diffusion along the two close-packed 
stepson Al (111). 



along the {100}-faceted step. Indeed, the calculations yield a value of 0.32 
eV, smaller than the value of 0.48 eV obtained for the hopping along the 
{lll}-faceted step. The calculated barrier for the exchange mechanism is 
about the same for both steps (0.42 and 0.44 eV). 

4.1.2. Ab initio KMC study of growth 

We now analyze typical growth conditions where kinetic processes are 
dominant. The detailed characterization of the energetics of diffusion pro- 
cesses carried out by Stumpf and Scheffler [|90| , riT^f| for A1/A1(111) and 
presented in the previous Section has provided several parameters for re- 
alistic KMC simulations. Among the processes listed in Table [T] we have 
considered the following diffusion mechanisms: 

(i) diffusion of a single adatom on the flat surface: E& = 0.04 eV; 

(ii) exchange diffusion from upper to lower terraces: E& = 0.06 eV at the 
{100}-faceted step and E& = 0.08 eV at the {lll}-faceted step; 

(in) diffusion parallel to the {100}-faceted step via hopping: E& = 0.32 
eV; 



[iv 



diffusion parallel to the {lll}-faceted step via exchange: E& = 0.42 
eV. 



The DFT calculations give that the binding energy of a dimer is 0.58 
eV [P0|1 , an d we therefore assume that dimers, once they are formed, are 



stable (i* = 1). Moreover, in the lack of reliable information we assume 
that dimers are immobile. We note that the reported value for the self- 
diffusion energy barrier is rather low (0.04 eV) |90| , |115|1 and comparable 
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to the energy of optical phonons of Al (111) (0.03 — 0.04 eV [|ll(j|1 ). Thus, 



simulations at room temperature may not be reliable because the concept 
of single jumps between nearest neighbor sites is no more valid. A single 
optical phonon can furnish enough energy to an adatom for leaving its 
adsorption site and diffusing on the flat surface. At room temperature the 
level population of optical phonon is high and the adatoms have practically 
no saddle point and migrate freely on the flat surface. We therefore limited 
our study to substrate temperatures T < 250 K. 

We adopt periodic boundary conditions, and our rectangular simulation 
area is compatible with the geometry of an fee (111) surface. The dimen- 
sions of the simulation area are 1718 x 2976 A 2 . These dimensions are a 
critical parameter and it is important to ensure that the simulation area 
is large enough so that artificial correlations of neighboring cells do not 
affect the growth patterns. The mean free path A of a diffusing adatom 
before it meets another adatom with possible formation of a nucleation 
center or is captured by existing islands should be much smaller than the 
linear dimension of the simulation cell [|117J1 . Since A is proportional to 
(D/F) 1 / 6 0, we find that (with F = 0.08 ML/s) A - 50 A for T = 50 



and gets as large as ~ 10 3 A for T = 250 K. We see that our cell is large 
enough (for the imposed deposition rate) for T < 150i^, whereas at higher 
temperatures the dimensions of the cell are too small, i.e., for T > 150i^ 
the island density is influenced by the sizes of the simulation area. Never- 
theless, the island shape is determined by local processes (edge diffusions) 
and is still meaningful. 

In the KMC program two additional insights extracted from the DFT 
calculations are included: (i) the attractive interaction between steps and 
single adatoms, and (ii) the fact that diffusion processes take place via 
different mechanisms (hopping or exchange). Particularly the second point 
plays an important role in our investigation. In several KMC simulations 
of epitaxial growth the attempt frequency of the diffusion rate has the 
same value for all the processes, and this value lies usually in the range 
of a typical optical phonon vibration or the Debye frequency. However, 
this assumption may not be right. First, processes with larger activation 
barriers may have a larger attempt frequency than processes with smaller 
energy barriers. This is a consequence of the compensation effect briefly 
described in Section |2.1| . Moreover, processes as hopping and exchange that 
involve a different number of particles and different bonding configurations 
may also be characterized by different attempt frequencies. This has been 
observed |8|, |118| - |rZ0|1 for several systems (Rh, Ir, Pt) and implies that the 
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Figure 14. A surface area of (1718 x 1488) A 2 (half of the simulation area) at four 
different substrate temperatures. The deposition rate was 0.08 ML/s and the coverage in 
each picture is 6 = 0.04 ML. 



attempt frequency for exchange diffusion can be larger by up to two orders 
of magnitude than that for hopping. For Al surfaces calculations with the 
embedded atom method Ijl^ljl showed a difference of prefactors of one order 
of magnitude. 

The results of the ab initio KMC simulations shown in Fig. [TJ] are for 
a coverage of G = 0.04 ML. When the growth temperature is 50 K the 
shape of the islands is highly irregular and indeed fractal. Adatoms which 
reach a step cannot leave it anymore and they cannot even diffuse along 
the step. Thus, at this temperature ramification takes place into random 
directions, and island formation can be understood in terms of the hit and 
stick model [T7j. At T = 150 K the island shapes are triangular with their 
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sides being {100}-faceted steps. Increasing the temperature to T — 200 K 
a transition from triangular to hexagonal shape occurs and for T = 250 K 
the islands become triangular again. However, at this temperature they 
are mainly bounded by {lll}-faceted steps. 

To understand the island shapes in the temperature regime between 150 
and 250 K we consider the mobility of the adatoms along the steps (at 
such temperatures the adatoms at the step edges cannot leave the steps): 
The lower the migration probability along a given step edge, the higher 
is the step roughness and the faster is the speed of advancement of this 
step edge. As a consequence, this step edge shortens and eventually it may 
even disappears. Since diffusion along the densely packed steps on the 
(111) surface (the {100} and {111} facets) is faster than along steps with 
any other orientation this criterion explains the presence of islands which 
are mainly bounded by {100}- or {lll}-faceted steps. The same argument 
can be extended to the diffusion along the two close-packed steps and 
applied to the triangular islands at T — 150 K, where the energy barrier 
for the diffusion along the {111} facet is larger and thus the {100}-faceted 
steps survive so that triangular islands with {100} sides are obtained. By 
considering the energy barriers we would expect only these islands, until 
the temperature regime for the thermal equilibrium is reached. However, as 
noted in Section [27]], the diffusion rates of adatoms are not only governed 
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Figure 15. Temperature dependence of the edge diffusion rates for atom diffusion along 
the {100}-faceted step by hopping with T = 2.5 x 10 12 s _1 (solid line), and along the 
{lll}-faceted step by exchange with T = 2.5 x 10 14 s _1 (dash-dotted line). 
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Figure 16. Shape of the islands at T = 250 K as they develop with time (or coverage). 
The snapshots refer to 6 = 0.006 ML, 6 = 0.02 ML, 6 = 0.025 ML, and 9 = 0.03 ML. 
The section of the simulation cell that is shown is 1718 x 1488 A 2 and the deposition rate 
is 0.08 ML/s. 



by the energy barrier but also by the effective attempt frequency. For 
A1/A1(111) the effective attempt frequencies have not been calculated, 



but the analysis of Ref. |5Uj proposes that the exchange process should 
have a larger attempt frequency than the hopping process. The results 



displayed in Fig. [P| are obtained with 1.0 x 10 s for the diffusion on 
the flat surface, 2.5 x 10 12 s^ 1 for the jump along the {100}-faceted step, and 
2.5 x 10 14 s" 1 for the exchange along the {lll}-faceted step. These effective 
attempt frequencies are the only input of the KMC not calculated explicitly 
by DFT, but were estimated from the theoretical PES as well as from 
experimental data for other systems. In Fig. [T5] the edge diffusion rates 



40 

along the two steps are plotted as a function of the reciprocal temperature. 
At lower temperatures the energy barrier dominates the diffusion rate but 
at T — 250 K the attempt frequencies start to play a role and lead to 
faster diffusion along the {111} facet than along the {100} one. Thus, 
the latter steps disappear and only triangles with {lll}-faceted sides are 
present. The roughly hexagonally shaped islands at T — 200 K are a 
consequence of the equal advancement speed for the two steps at that 
temperature. Obviously, the temperature dependence of the growth shapes 
found in Fig. [P| is crucially determined by the ratio of the two diffusivities 
and in particular by the temperature at which the two lines of Fig. |T5| 
cross. If the difference were only one order of magnitude, the crossing 
would be at a temperature that is too high (500 K). The formation of 
fractals (Fig. [T3], upper left) and of {100}-faceted step triangles would 
still occur. Obviously, the importance of the attempt frequencies should 
receive a better assessment through accurate calculations, and work in this 
direction is in progress. 
A peculiarity of the triangular islands in Fig. \TQ is that they exhibit con- 



cave sides. In order to understand this behavior we examine the evolution 
of the island shape for the deposition at T = 250 K. The results are col- 
lected in Fig. [16]. At very low coverage the islands are roughly hexagonal 
and upon successive deposition they evolve into a nearly triangular shape. 
The longer sides are formed by straight {lll}-faceted step edges but short 
{100}-faceted edges can still be identified, at least for O < 0.01 ML. The 
latter edges become rougher and progressively disappear. For O = 0.025 
ML the sides are still nearly straight, but at O = 0.03 ML the concav- 
ities appear. The corners of the triangles seem to increase their rate of 
advancement during deposition. The effect can be understood on the ba- 
sis of competition between adatom supply from the flat surface and mass 
transport along the sides. The adatom concentration field around an island 
exhibits the steepest gradient close to the corners, and the corners of the 
islands receive an increased flux of adatoms. When the sides of the islands 
are not too long, this additional supply of adatoms is compensated by the 
mass transport along the steps, i.e., the adatoms have a high probability 
to leave the region around the corners before the arrival of the successive 
adatom. For O = 0.025 ML this scenario still seems to be true, while at 
O = 0.03 ML the island edges are longer and the mass transport along the 
sides is not able to compensate the additional supply of particles at the 
corners. That means that the probability for a particle to leave the corner 
region and to move along the island edge before being reached by another 
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particle decreases considerably, and the corners start to grow faster than 
the sides of the triangles so that the concave shape develops. 

4.2. Ag(lll) 

4.2.1. The influence of strain on surface diffusion 

Growth of one material on a different material is of particular interest for 
a number of technological applications. In such a heteroepitaxial system 
with usually different lattice constants the material to be deposited is under 
the influence of epitaxial strain. Growth of Ag on Pt (111) and Ag on a 
thin Ag film on Pt (111) has been the focus of a number of recent studies 
[|112| , |122| , |123|1 , and with a lattice mismatch of 4.2% it serves as an ideal 
system that can provide important information about the effects of strain 
during growth. We will particularly discuss how strain affects the surface 
diffusion barrier. 

Only few theoretical studies of the effect of lattice mismatch on the diffu- 
sion barrier are present in the literature. For a metallic system we are only 
aware of results for Ag on Ag (111) where the authors of Ref. [|112|1 find in 



an EMT calculation that the diffusion barrier increases under tensile strain 
and decreases under compressive strain. 

Here, we present first principle calculations (more details are given in 
Ref. fl52fl ) where we study systematically the dependence of the diffusion 



barrier on the lattice constant for Ag on Ag(lll) [|112| j. In the range of 
±5 % strain the DFT results exhibit a linear dependence with a slope of 
0.7 eV as it is illustrated in Fig. [17]. The calculated diffusion barrier for 
the unstrained system, E d g ~ g = 81 meV, is in good agreement (within the 
error margins of the experiment and the calculations) with the scanning 
tunneling microscopy (STM) results of E^ g ~ Ag = 97 meV. The accordance 
between experiment and theory extends to the system Ag/Pt(lll) and 
Ag/IML Ag/Pt (111). These results are summarized in Table g. In Fig. [T7j 
the DFT-LDA results are compared to those of an EMT study [|TT211 . The 



EMT results exhibit a linear dependence only for very small values of strain 
(±2 %) and the diffusion barrier starts to decrease for values of misfit larger 
than 3 %. Indeed, it is plausible that a decrease of the diffusion barrier 
occurs when the atoms are separated far enough that eventually bonds 
are broken. However, as our DFT-LDA results show, for Ag/Ag (111) this 
happens at values for the misfit that are larger than 5 %. Additionally, 
when comparison with experiment is possible [i.e., Ag on Ag (HI), and Ag 
on a monolayer Ag on Pt (HI)] the EMT results are off by a factor that 
varies from 1.2 to 2. 
The DFT results in Fig. |T7| were obtained with the LDA for the exchange- 



42 



120 



20 L- 1 



~i ■ ■ ' 



T 1 1 T" 




^ 



0.95 1.00 

Relative Lattice Constant 



1.05 



Figure 17. Diffusion barrier (in meV) for Ag on Ag(lll) as a function of strain. The 
circle s are DFT-LDA results from Ref. |)2|] and the squares are EMT results from Ref. 



correlation functional and test calculations show that GGA increases the 
diffusion barrier by no more than 5 — 10%. The (lll)-surface is a closed 
packed surface with a very small surface corrugation and since LDA and 
GGA results on this surface do not show significant differences for the dif- 
fusion barrier (as long as the mechanism is hopping and not exchange) it 
is plausible to assume that both the LDA and the GGA are good approx- 
imations for the exact exchange-correlation functional. This is also true 



for Pt/Pt(lll) [PI and Ag on Ag(100) g?|. The general trend of an 



increasing energy barrier for hopping diffusion with increasing lattice con- 
stant is quite plausible (for exchange diffusion see Section |4.4.2|) . Smaller 
lattice constants correspond to a reduced corrugation of the surface, and as 
result the atom is not bonded much stronger at the adsorption sites than at 
the bridge site. In contrast, when the surface is stretched the corrugation 



Table 2 

Diffusion barriers (in meV) for Ag on Pt (111), Ag on one monolayer (ML) Ag on Pt (111) 
and Ag on Ag(lll). 



System 



Ag/pt(in; 

Ag/1ML Ag/Pt(lll) 

Ag/Ag(lll) 



Experiment ||TT2|| EMT ||TT2|| DFT §2 



157 

60 

97 



81 
50 
67 



150 
65 

81 
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increases and the adsorption energy at the three-fold coordinated hollow 
sites increases. This picture will change when the strain is so large that 
bonds are broken and then it is expected that the hopping diffusion barrier 
will start to decrease again at very large tensile strain. 

It is worth noting that the diffusion barrier for Ag on top of a pseudo- 
morphic layer of Ag on Pt (111) is substantially lower than it is for Ag on 
Ag(lll). A question that arises is whether this reduced diffusion barrier 
is a result of the compressive strain or should be ascribed to electronic 
rearrangements induced by the Pt substrate. The diffusion barrier for 
Ag on Ag(lll) with a lattice constant that is compressed to the value 
of the lattice constant for Pt is E^ g ' Ag = 60 meV while that for Ag on 
Pt (111) (also with the Pt lattice constant of 3.92 A obtained from DFT) 
is E^ g ' Ag/Pt — 65 meV. The agreement of these two values suggests that the 
reduction of the diffusion barrier for Ag on a layer of Ag on Pt (111) is 
mainly a strain effect and that the diffusion barrier on top of a layer of Ag 
is essentially independent of the substrate underneath. 

Brune et al. [ |112j| also measured the island densities of Ag on two mono- 
layers (ML) of Ag on Pt (111) and found that the island density is much 
larger than it is for Ag on just one ML of Ag on Pt (111). But the reason 
for this increased island density is not a larger barrier for surface diffusion. 
The second layer of Ag on Pt (111) reconstructs in a trigonal network where 
domains with atoms in the fee and hep site alternate JT25J. This recon- 



struction occurs either during growth with high enough adatom mobility 
or upon annealing and it can be concluded that this trigonal network is the 
equilibrium structure. The periodicity of these domains is approximately 
two domain boundaries for every 24 atoms. This can be understood very 
well with purely geometrical arguments because the lattice mismatch is 
4.2 % and every domain boundary implies that there is half of an Ag atom 
less so that the domain network provides an efficient mechanism to relief 
epitaxial strain. The barrier to diffuse across such a domain wall appears 
to be rather high and domain walls act as repulsive walls so that the island 
density is determined by the defect density and not the barrier for self 
diffusion. It is not clear however why this domain network is formed only 
after 2 ML Ag have been deposited and not already upon completion of 
the first Ag layer. 

To answer this question the adsorption energy of an adatom in the fee 
and in the hep site were compared by Ratsch et al. [|52fl . Calculations were 
carried out with a (1 x 1) and a (2 x 2) cell and slab thicknesses of up to 
5 layers. It was found that the fee site is energetically more favorable than 
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Table 3 

Energy difference AE = —(Ef cc — -Ehcp) between the total energies of the two adsorption 
sites for Ag on Pt (111) and Ag on Ag (111). 

System AE (in meV) AE (in meV) 

__^__ 1 x 1 cell 2x2 cell 

Ag/Pt(lll) 30 50 

Ag/IML Ag/Pt(lll) - <5 

Ag/Ag(lll) <10 <10 



the hep site in all cases. The energy differences between the two adsorption 
sites for the first and second Ag layer on top of Pt (111) and Ag on Ag (111) 
are summarized in Table |3|. During deposition of the first layer the fee site 
is energetically preferred compared to the hep site, and the film growth 
pseudomorphically with the atoms in the fee sites. However, after the first 
layer has been completed, the total energies of these two sites are almost 
indistinguishable so that the system is not prevented from reconstructing 
in the described domain network to relieve epitaxial strain. 

4.2.2. The role of antimony as a surfactant 

As already mentioned, Ag atoms deposited on Ag(lll) form "mountains" 
(or "mounds" ) as it has been seen by RHEED |126| , x-ray reflectivity ex- 



periments [|114j| , and STM [|g^Jg3|JI27| for a wide range of temperatures. 



The situation changes completely when small amounts of antimony (0.2 
ML) are deposited on the surface and the growth mode becomes two di- 
mensional [|4T3|,|T27|1. 

The natural question that arises is: What is the action of Sb? These 
contaminants are called surfactants, although this term may be misleading. 
In the original definition a surfactant should reduce the surface energy. 
Particularly for metallic systems, however, the contaminants rather affect 
the kinetics of the processes and change the growth mode. 

In Section \Z.'J\ we have described some possible scenarios for the action 



of these surfactants. Clearly, different surfactants work by different mech- 
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Figure 18. Adsorption sites of Sb on Ag (111) considered in the calculations of Ref. ||128 |. 



anisms (or a combination of them). However, as already pointed out in 
Section |2.3| a good surfactant has to satisfy one essential requirement: It 
has to stay on the surface. To investigate this crucial point, Oppo et 
al. pn~1 , |128j| carried out DFT-LDA calculations for the adsorption energies 
of Sb atoms located in different positions on Ag(lll). These calculations 
have given important insight into the action of Sb during growth. The 
evaluated binding energies for Sb are collected in Table f| for different cov- 
erages. It is clear that the substitutional site is greatly favored with respect 
to on-surface fee adsorption and sublayer adsorption for all coverages con- 
sidered (Osb = 1/4, 1/3, and 1 were considered)). Sb is thus expected to 
be adsorbed in substitutional sites at not too high Sb coverages and this 
geometry is shown in Fig. [TR An incorporation of the contaminant into 



the Ag bulk can be ruled out because its energy cost is too high. The Sb 
atoms segregate to the surface and do not incorporate into the bulk, since 
its sizes are somewhat too large for a bulk vacancy, but appropriate to fit 
into a surface vacancy. 

The DFT results clarify the energetical ranking for the adsorbed Sb atom, 
but the understanding of its action on the Ag atoms requires the determi- 
nation of its influence on the Ag adsorption. Two possible locations of the 
Ag adatom on the substitutional Sb-covered surface are considered: a near 
site and a far site, depending on whether Ag and Sb are nearest neighbors 
or not. These two sites are labeled 3 and 2 in Fig. [T8|, and their calculated 
adsorption energies are 1.99 eV/atom and 2.02 eV/atom, respectively. The 
Ag adatom prefers energetically to sit on the clean portion of the surface 
(clean site with i? ad g s = 2.41 eV/atom), whereas close to an incorporated 
Sb atom it favors the far site. Thus, Ag adatoms avoid the vicinity of the 
Sb atoms. 
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A schematic summary of the action of Sb during growth of Ag on Ag (111) 
is given in Fig. [TO|. The substitutional antimony atoms act as repulsive cen- 
ters for the silver adatoms (upper panel), whose mean free path is reduced 
by the presence of the contaminants (middle panel). Thus, small-sized Ag 
islands appear on the surface with an increased density. Moreover, the 
antimony atoms are expected to affect the form of the Ag islands, and 
the islands will grow in directions where they can avoid to get close to 
Sb centers. Finally, just before the Ag-Sb surface alloy layer is covered 
by the newly deposited Ag layer, the Sb atoms find themselves in the 
disfavored subsurface configuration (see, for instance, Table [|), and seg- 
regate to the new surface layer (Fig. |T9|, lower panel). The alloy layer is 
now reestablished and the process can start again. The combination of all 
these features yields layer-by-layer growth. It is interesting to note that 
according to Eq. (||) only an increase of island density by a factor of 2 is 
necessary to induce two dimensional growth. Vrijmoeth et al. [|127|] report a 
much higher increase of the island density (a factor of ~ 6) of the annealed 
surface in the presence of Sb. 

A recent theoretical study JT29J where a rate-equation analysis has been 





Figure 19. Sketch of the action of Sb on Ag(lll) 
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combined with Monte Carlo simulations has confirmed the main idea pro- 
posed in the work of Oppo et al. [ gT| , [I2gfl . The main role of the Sb repulsive 
network is to reduce the mean free path of the Ag adatoms, and as a re- 
sult the island density is enhanced. The changed character of diffusion 
leads to new scaling relations in which the island density shows a strong 
dependence on the impurity concentration. 

4.3. Microscopic processes at Al (100) 

Self-diffusion at flat regions of close-packed surfaces is typically consid- 
ered to proceed by series of hops between adjacent equilibrium adsorption 
sites. For a fee (100) surface an adatom in a fourfold hollow site moves over 
the twofold bridge site to a neighboring fourfold hollow site [see Fig. |5](a)]. 
After some time an adatom has visited most fourfold hollow surface sites, 
and the visited sites form a (1 x 1) pattern. An alternative mechanism 
for surface diffusion is atomic exchange where the adatom moves by dis- 
placing a neighboring surface atom. On a fee (100) surface the exchange 
process occurs along the [010] and [001] directions. The transition state 
[see Fig. |5|(b)] may be described as a dimer, consisting of the adatom and 
the lift-up surface atom located above a surface vacancy. Subsequently the 
displaced surface atom becomes a new adatom in a next-nearest-neighbor 
fourfold site. When this diffusion mechanism is active, an adatom visits 
only every second fourfold hollow site at the surface, and the visited sites 
form a c(2 x 2) pattern. This difference in the patterns of visited sites 
represents clear evidence of the active mechanism. Indeed, in this way self- 
diffusion by exchange has been experimentally observed by Kellogg and 
Feibelman for Pt (100) [|TTJ and by Chen and Tsong for Ir (100) |jI2j|. 

A theoretical study of exchange diffusion was performed by Feibelman [jT^] 
for Al (100). Feibelman performed DFT calculations and predicted a rather 
high energy barrier for hopping (E d op = 0.65 eV) in conflict with exper- 
imental estimates. Inspired by ideas of Pandey [|25f Feibelman realized 



that at surfaces it is more appropriate to think of self- diffusion in terms 
of making and breaking of chemical bonds rather than in terms of a hard 
sphere rolling over a bumpy plane. In fact, his calculations of the exchange 
mechanism show that its barrier is substantially smaller (by ~ 0.45 eV) 
than that of the hopping process. Thus, he predicted that for Al(100) 
self-diffusion proceeds by atomic exchange. The process is caused by the 
noticeable covalency of aluminum, which can form directional bonds (by 
sp hybridization) at certain atomic geometries. In the transition state of 
exchange diffusion at Al(100) [see Fig. 0(b)] each atom of the dimer that 
sits over the vacant site forms three chemical bonds, two with surface Al 
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atoms and one with the other atom of the dimer. Inspection of the elec- 
tron density |T3j revealed the formation of directional bonds and that the 



coordination of the surface atoms at the transition state can be described 
as three-fold. In other words: At the transition state of hopping diffusion 
the adatom is two-fold coordinated, while at the transition state of the ex- 
change diffusion it is three-fold coordinated. For a group III atom it thus 
appears plausible that the latter transition state has a lower energy. We 
note in passing that recent embedded atom calculations have questioned 
the importance of the exchange mechanism at Al (100) [|l2l|| , but DFT cal- 



culations of Stumpf and Scheffler [J9D| fully confirmed the scenario proposed 
by Feibelman. 

These studies are predictions with only indirect experimental support, 
namely the observation that the diffusion barrier is low. Only for two 
other fee (100) surfaces exchange diffusion has been identified [Pt (100) and 
Ir(100)]. Apparently it does not occur at other fee (100) transition-metal 
surfaces. Recent calculations have shown that the mechanism that stabi- 
lized the exchange diffusion at the late hd transition metals is significantly 
different from the covalency mechanism at Al(100). In Section |4.4.^| we 
will summarize these results. 

4.4. Ag(100) 

4.4.1. Microscopic processes 

For silver it is known that growth on the (111) surface is three dimen- 
sional because a noticeable step-edge barrier exists that has been estimated 
by various authors using STM |3^,[33j to be between 0.1 and 0.15 eV. This 
additional step-edge barrier explains why the growth of silver perpendicu- 
lar to the (111) surface proceeds in a three-dimensional mode. 

In order to analyze why the growth of the Ag (100) surface is qualita- 
tively different (two dimensional) Yu and Scheffler [gBj performed DFT 
calculations within the LDA and the GGA, and the results are listed in 
Table 0. In contrast to the above discussion for Al(100) self-diffusion on 
Ag(100) proceeds by hopping. Exchange diffusion has an energy barrier 
which is about 0.3 eV higher. 

However, according to the DFT calculations diffusion across a descending 
close-packed step takes place via the exchange mechanism (cf. Table [5]). 
The calculated activation energy for the step-down motion by exchange is 
(within the numerical accuracy) identical to that of the hopping diffusion 
at the terrace. Figure pO shows the total energy curves along the diffusion 



paths for hopping and exchange diffusion. 
The low energy barrier of the exchange diffusion can be understood by 
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considering the geometries along the reaction pathway: Atom 2 moves out 
of the step edge and the process is energetically stabilized because adatom 

1 keeps in contact with atom 2 and finally occupies the original site of atom 

2 (Fig. ^D], right side). The transition state for the exchange process is near 
the bridge site formed by two step-bottom atoms on the lower terrace (S ex 
in Fig. £0|, right side). Inspection of the geometry at S ex (see Ref. |86j for 
details) shows that the local coordination of atoms 1 and 2 remains high. 
At S ex each of theses atoms is five-fold coordinated. This finding that 
there is no additional energy barrier provides a natural explanation for the 
smooth 2-D growth of Ag (100). By analyzing the STM images of Ag(100) 
during deposition Zhang et al. [|l3L)j| gave an estimate of 0.025±0.005 eV for 
the additional step-edge barrier in good agreement with the DFT results. 
We also note that the coordination number arguments employed above are 
valid in general, thus we expect that our finding that step-down diffusion at 
Ag(100) proceeds by the exchange mechanism also applies to other noble 
and other fee (100) transition-metal surfaces. 

Yu and Scheffler also studied the diffusion parallel to step edges and 



the step formation energies [M]. They find that diffusion parallel to step 



edges is faster than that at flat surface regions and that the closest-packed 
step, the (110)/{111} step dominates. Thus, it is predicted that step 
edges are rather straight and that in equilibrium islands have an approx- 
imately square shape. To be precise, the size is that of an octagon with 
long (110)/{111} and shorter (100)/{110} steps. From the step formation 
energy the length ratio was determined to be 10:3. 

4.4.2. The influence of strain on surface diffusion 

As noted above, self-diffusion on fee (100) surfaces may proceed by the 
hopping process (as for Ag) or by exchange (Pt, Ir, and probably also Al). 
Yu and Scheffler [jnj] recently predicted that exchange diffusion should 
occur on Au(100) and strained Ag(100), and two mechanisms were em- 



Table 5 

Energy barriers E& (in eV) for different processes on Ag (100) calculated within the LDA 
(GGA). 

process mechanism E& (eV) 

flat Ag(100) hopping 0.52 (0.45) 

flat Ag(100) exchange 0.93 (0.73) 

step || hopping 0.30 (0.27) 

step _L descent. exchange 0.52 (0.45) 

step _l_ descent. hopping 0.70 (0.55) 
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phasized that stabilize the exchange over the hopping mechanism: i) the 
tensile surface stress, and ii) the correlation of bond strength and local 
coordination. 
For Au(100) the rather low energy of the exchange-diffusion transition 



state 



[Ef == 0.65 eV, 



ihop 



E d p = 0.83 eV) could be correlated with its geometry 
(compare Fig. 0(b)): The dimer is only 1.29 A above the surface layer, 
which is 37 % closer to the center of the top layer than the inter-layer 
spacing in the bulk. Thus, the two atoms of the dimer interact with the 
atoms of the top layer of the substrate but get rather close to the second 
layer. In other words, the surface likes to attain closer packing which 
reflects its noticeable tensile stress [ |T5Tj| . Indeed, the surface stress of 
the late 5d metals is significantly higher than that of their 4d and 3d 
isoelectronic elements. The difference has been traced back to relativistic 



effects | |lal|1 , which play a noticeable role for the heavier 5d metals: The 



relativistic effects give rise to a contraction and energy lowering of s-states 
and as a consequence, the d-band moves closer to a Fermi energy [|l3l 



Indeed, a relativistic treatment is most important in order to attain a 
good description of structural and elastic properties of 5d metals, while 
it is not important for the 4d metals. Thus, while the significant tensile 




OO 

OOt 

OO 1 ' 




OO 
OO 

oo, 

OO 1 ' 



„ 1.0 

> 

CD 



CD 
LU 



0.5 



0.0 




0.0 



1.0 2.0 3.0 

Dis1ance( a s ) 



4.0 



10 



> 

CD 



£l 0.5 

CD 

w 00 



- 






- 






- 


\J 


\ 


I 


I 


V 

M 

i , i 



-1.0 0.0 1.0 2.0 

Dis1ance( a s ) 



Figure 20. Diffusion path (upper panel) and total energy (lower panel) of an Ag 
adatom diffusing across a descending step by hopping (left side) and exchange process 
(right side). The values of the energy have been calculated within the LDA. For the 
exchange mechanism the total energy is plotted as function of the distance of the step 
edge atom 2 from the undistorted step edge. 
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surface stress of Au (100) pulls the dimer of the exchange transition state 
"into" the surface, it lowers the energy of the transition state, and enables 
exchange diffusion. The surface stress at Ag(100) is too weak to have a 
significant effect. 

With respect to the stress a stretched silver film gets more gold-like, and 
therefore the two diffusion mechanisms for a strained silver slab were ana- 
lyzed in detail. When the system is under tensile strain the energy barrier 
for the exchange diffusion increases, while the barrier for the hopping dif- 
fusion decreases. For hopping diffusion the trend can be understood as fol- 
lows: Smaller lattice constants correspond to a reduced corrugation of the 
surface potential, and thus diffusion energy barriers are reduced |62| , |ll2J| . 
In contrast, for a stretched surface the corrugation increases and the ad- 
sorption energy at the four-fold coordinated hollow sites increases. The 
latter reflects the wish to reduce (at least locally) the strain induced sur- 
face stress. The hopping-diffusion transition state is less affected by the 
strain than the adsorption site. For exchange diffusion this is just the op- 
posite. Here the transition-state geometry reacts particularly strongly to 
the tensile stress, and locally the tensile stress is reduced by the very close 
approach of the dimer [Fig. 0(b)] . Thus, it is predicted that for pseudo- 
morphic Ag films (with increased parallel lattice constant) self-diffusion 
should get noticeably affected by the exchange mechanism. The results of 
this study strongly suggest that tensile surface stress (to be precize, the 
excess surface stress, i.e. the strain derivative of the surface energy) is the 
main actuator for the exchange diffusion on fee (100) surfaces [[T4"j| . 
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